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Abstract 

This paper establishes a global contraction property for networks of phase-coupled oscillators characterized 
by a monotone coupling function. The contraction measure is a total variation distance. The contraction 
property determines the asymptotic behavior of the network, which is either finite-time synchronization or 
asymptotic convergence to a splay state. 
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1. Introduction 

Networks of coupled oscillators are a general paradigm to understand a wealth of natural phenomena 
@> O EH as well as an efficient model for the design of engineered systems 0, EH • In the limit of weak 
coupling, realistic models of limit-cycle oscillators (evolving in high-dimensional spaces) can be reduced to 
one-dimensional models of phase oscillators (evolving on the circle) coupled through their phase differences 
Though more amenable to mathematical analysis, these generic phase-coupled models may exhibit 
rich and complex ensemble behaviors and have attracted intense research interest in the past decades (e.g. 
Kuramoto model [||). 

A phase-coupled model is characterized by its coupling function (which is closely related to the phase 
response curve of the oscillators) . The seminal work of Kuramoto [B[ assumes a sinusoidal coupling function. 
The present paper rather considers a network of oscillators characterized by a monotone coupling function 
(we denote these oscillators as monotone oscillators). This model was first studied in [fjj], through the phase 
reduction of the popular leaky integrate-and-fire neuron model. Using local analysis, the author investigated 
the effect of a delay on the network stability and showed that the periodic collective motion of the oscillators 
becomes quasiperiodic when the delay exceeds a critical value. Since it is known that delayed oscillators 
lose their monotonicity property, this result suggests a link between the stability properties of the oscillators 
and their monotonicity properties. 

In contrast to the local results provided in Q , we present in this paper a global contraction property for 
networks of monotone oscillators. The result is shown with respect to a 1-norm which has the interpretation 
of a total variation distance and which is inspired from our previous studies 8, 9]. The contraction property 
of the model determines the asymptotic behavior of the network. Namely, monotone oscillators either achieve 
perfect synchronization in finite time or converge to a unique anti-synchronized state (splay configuration). 

The paper is organized as follows. Section [5] introduces the model of monotone phase-coupled oscillators. 
In Section^ we present our main result on the contraction property of the oscillators. The result is exploited 
in Section 5] to study the collective behaviors of the network. Finally, the paper closes with some concluding 
remarks in Section [5] 



2. Monotone phase-coupled oscillators 

We consider a network of identical phase-coupled oscillators. A single (uncoupled) phase oscillator k is 
characterized by a phase 9k £ S 1 (0, 2ir) that evolves on the circle with constant velocity 0\ = cu, where ui 
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is the natural frequency of the oscillator. Within the network, N phase oscillators are (all-to-all) coupled 
through their phase differences: they evolve on the TV-torus T N = S 1 x • • ■ x S 1 according to the canonical 
phase dynamics 

N 

9 k =u + Y,r(9k-6i) k = l,...,N. (1) 

j'=i 

Using averaging techniques, every network of weakly-coupled identical limit-cycle oscillators can be re- 
duced to the form ((T]), with an appropriate coupling function T(-) that is closely related to the phase response 
curve of the oscillators 0,0]. The phase response curve characterizes the phase sensitivity of the oscillator to 
a (infinitesimal) perturbation, a quantity that can be either numerically computed or experimentally mea- 
sured. Since it is a function computed on the (periodic) limit cycle, the phase response curve is 27r-periodic, 
and so is the coupling function T(-). 

We assume that the coupling function satisfies the following monotonicity assumption, in which case the 
phase-coupled oscillators are called "monotone oscillators". 

Assumption 1 (Monotonicity). The coupling function T(-) is strictly monotone, i.e. either T'{9) > 
V0 E (0,2tt) or r'(0) < ^9 e (0,2tt), with r denoting the first derivative ofT. 

Remark 2.1. Assumption Q] implies that the coupling function is discontinuous, that is r(0~) ^ r(0 + ). 
However, it is important to note that the 27r-periodicity condition imposes T(0) = r(27r) (Figure [TJ. 




The monotonicity property is frequently encountered in spiking oscillators, at least in good approxima- 
tion. Monotone phase-coupled oscillators were first obtained in the case of (weakly) pulse-coupled leaky 
integrate-and-fire oscillators Q. More generally, for a purely impulsive coupling, monotone oscillators ap- 
pear as the phase reduction of oscillators characterized by a monotone phase response curve, such as van 
der Pol oscillators with strong relaxation or limit-cycle oscillators near a homoclinic bifurcation [l[ . 



Remark 2.2. The analysis of a network of monotone oscillators is well-known when all oscillators are 
initialized within a semicircle: in this case, a change of coordinates maps the dynamics into M. N and the 
model falls into the well-studied category of consensus models on a convex set, see e.g. [ll|, 14 1. In the 
present paper, the emphasis is on the whole (non-convex) iV-torus T N . 



3. Main result 

Definition 1 (Contraction). Let S be a continuous-time dynamical system defined on a metric space X 
(with the distance d) and let $ denote the flow associated with S, i.e. 4>(x, t) is an orbit of S, with 

2 



4>(x, 0) = x G X. Then S is contracting (resp. expanding) in Xq C X with respect to if 

^d(#(x,£),*(y,£)) <0 (resp. >0) 

for all x ^ y and for all £ such that 4>(x, £) G Xq, 3?(y, £) elo- 

We will study the contraction of model (Q} with respect to the distance 

n-1 

d(x,y) = ||x - y|| (1) = |xi -yi\ + ^2 \xj+i - y j+ i - (xj - %)| + |x n - y n \ . (2) 

j'=i 

This particular 1-norm distance is motivated by previous results on intcgratc-and-firc models. For finite 
populations of leaky integrate-and-fire oscillators, the distance ([2J was successfully used in Q to capture 
the contraction property of the so-called firing map (a discrete-time map that provides snapshots of the 
network configuration at the successive firings of the oscillators [10]). For infinite populations, the continuous 
equivalent of ^ induces a Lyapunov function for the partial derivative equation that governs the evolution 
of the population density 0]. 

The distance has the remarkable interpretation of a total variation distance. Indeed, ^ corresponds 
to the total variation of a piecewise linear function that interpolates the values xi — yi (Figure [5]) . In the 
continuous case (related to infinite populations) , this interpretation still holds since the continuous equivalent 
of is the Li-norm of the derivative 0|. 




Figure 2: The distance ((JJ is interpreted as the total variation of the linear interpolation of the values Xi — j/j. In the example 
of the figure, the distance is given by ||x — y||(i) = 2Ai + 2A2, a quantity which corresponds to the total variation of the 
piecewise linear function. 



Before establishing our main result, we must notice that the dynamics (JTJ) are invariant with respect 
to a rigid rotation of the oscillators, a property which prevents a contraction in the full space T^. We 
therefore remove this marginally stable rigid mode and consider the equivalent dynamics expressed in a 
rotating frame associated with an oscillator (without loss of generality, we choose oscillator 1). Denoting 
the phase differences by 9k = 6k+i — #i £ [0, 2ir] for k = 1, . . . , N — 1, the dynamics ((T|) arc rewritten as the 
(N — l)-dimensional dynamics 

N-l N-l 

Ok = r(e fc ) + r (^- - s i) - E r (-^') ■ ( 3 ) 
j=i 3=1 

Without loss of generality, we assume that the oscillators satisfy the phase ordering 9^ < 0k+i, which does 
not change over time since the oscillators are identical. Solutions ©(£) = (#i(£), ■ ■ ■ , 0N-i(t)) of © do not 
evolve in T^ 1 , but only in the closure C of the open cone C = {& G (0, 2?r) iV '- 1 |l9 fe < 9 k+1 , k = l...,N-2}. 
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An orbit that reaches the boundary of C corresponds to the synchronization of at least two oscillators. 

Now, under a mild technical assumption on the curvature of the coupling function, we are in position to 
prove the contraction property of the dynamics ([3]) with respect to the 1-norm distance The result is 
summarized in the following theorem. 

Theorem 3.1. Consider a coupling function that satisfies (i) Assumption^ (monotonicity) and (ii) either 
Y"{6) > V0 G (0,2tt) or T"{6) < V0 € (0,2tt). Then, © is either contracting (ifV < 0) or expanding 
(ifV>0)inC with respect to © . 

Proof. We prove that the dynamics arc contracting when T' < 0. Without loss of generality, we suppose 
that r" < and, considering the two orbits = (6x,-.., On-i) an d * = (ipi, . . . , ipN-i), we assume that 
#1 — i>\ < 0. The proof of the other cases follows on similar lines. 

We consider the distance © between and and remove the absolute values to obtain 




(4) 



k=l 



The subscripts K,(k) correspond to N c critical oscillators (1 < N c < N — 1) whose phases are characterized 
by a change of the sign of 6j+x — 4>j+i — (fij — V'j)- Formally, the map JC : {1, . . . , N c } n- {1, . . . , N — 1} is 
defined so that K.(k — 1) < JC(k) and so that 



[^/C(fe)+l - i>lC(k)+l - (^/C(fe) - ^JC(fc))] \PK(k) - i>K{k) - (0K(k)-X ~ ^K{k)-l)] < 0, 



with 9q — -00 = ®n — i>N = 0. The index /C(fc) corresponds to a maximum of 9j — tpj when k is even and to 
a minimum when k is odd (Figure [3]). In particular, one has 




K (k)-1 K(k) K (k)+1 j 




K (k)-1 K(k) K (k)+1 j 



(a) k even 



(b) k odd 



Figure 3: The critical oscillators correspond to extremal values of 9j — if)j . 



The time derivative of (HJ) yields 




(6) 



Next, we consider separately each term T^ k \ It follows from ([3]) that 



r « = (_i)* Y{§ K{h) )-T$ K{h) )+ r^w-^O-rfe)-^-) 



-^rc-^-rc-^-) 
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Using the mean value theorem, we obtain 



A rp{k) 
— 1 K{k) 

N-l 

+ E r 'fe( fc)J )(-i) fc (h(k) - Mk) - fa - fa)) 

._ 1 v v ' 

+ (-i)M-E r (-^)- r (-^)j' 



with 6c(fc) G [0jc(k)i or € fc(k)>0/C(k)] and £/C(k),j G [(0/c(k) - 0j) mod 2tt, ~ mod 27r ] or 

G [(^(fc) — V'j) mod 27r, (6ic(k) — Oj) mod 271"]. Since T" < 0, the values satisfy 

< Gc(k+i) fc < iv c , (8) 

C)C(fc),j < &e(k+i),j fc<iV c , j < K(k) or j > + 1) , (9) 

^K(N c ),j < &c(i)j IC(D < 3 < HNc) ■ (10) 

Because of the discontinuity of T, these relationships are not valid for an orbit that lies in the boundary of C. 

Disregarding the terms Ts, we show that each positive term of (in ||2J)) can be associated with 
a term T^ k+1) of T^ k+1 ) so that the addition of both terms T^ k) + T^ k+1) is negative. 



• If the term T^.A > 0, it is associated with the term T^ k+1 y It follows from (JSJ and from T" < that 



KYfc) > 0, it is associated with the term 
r'(6c(fc+i)) < r'(^(fe)) < 0. In addition, © implies 

-(-l) fc+1 (0jc(*+i) ~ fe(k+i)) < (-1)" (*jc(k) ~ ^(fc)) < ■ 
Then, we obtain ig^ + rg+ft, < 0. 

• If the term T^. > for any j 7^ K,(k), it is associated with the term +1 ' , j ^= )C(k + 1). (According 

to ([5]) and the definition of the critical oscillators, one has never rj fe) > for K(k) < j < K(k + 1).) 
It follows from © and from T" < that r'(^ (A . +1 ) J ) < r'(^(fc)j) < 0. In addition, © implies 

-(-i) k+1 (e K(k+1) - i, K{k+l) - fa - < (-i) k (e K(k) - i> K{k) - fa - v^O) < 0. 

Then, we obtain T^ k) + T^ k+1) < 0. 

For k — N c , the terms of cannot be associated with the terms of T^ k+1 \ Then, it remains to 
consider the terms T^ Nc ^ and the terms ( — l) k T^. We distinguish two cases: the case N c even and the case 
N c odd. 

Case N c even. The addition of the iV c terms (-l) fe T s yields EfciiM)^ = °- In addition, the 
relationship ([5]) implies 

^K(l) ~ 4>K(1) < 0K(N c ) ~ i>K(N c ) > (11) 

so that the term Tj^Jj^ > < 0. 
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If the term T^"' > for any j ^ K,(N C ), it is associated with the term , j ^ /C(l). (According to 

(fTTj) and the definition of the critical oscillators, one has never T^ Nc] > for j < £(1) or for j > IC{N C ).) 
It follows from (JTUJ) and from T" < that T'(^ic(x),j) < r'(£/e(jv c ) j) < 0. In addition, the inequalities (fTTj) 
imply 

-(-l) 1 fei) - ^(i) - & - 4'j)) < (-1)^ (*/cwo - &w - & - &)) < • 

Then, we obtain T^ Nc) + t] 1} < 0. 

Case N c odd. The addition of the iV c terms (— l) fe T s yields 



N c N-l 

5>i) fe r s = -t e = - £ r'K-i) $ - 
fe=i j=i 

with £_j G [27r — 27r — ipj] or G [27r — xpj, 2ir — 9j]. Since T" < 0, the values satisfy 

6c(iv )j < J < W) , (12) 

e-j < e^wj i > ^(i) ■ (13) 

In addition, the relationship ([5]) implies 

0JC(1) ~ i>K(i) < ^/C(JV C ) ™ ^K(N C ) < (14) 

so that the term T^f^ > < 0. 

• If the term T^- ^ > for any j 7^ K.(N C ), it is associated with the term —T'(j^-j)(6j — 'ipj)- (According 
to (fT4"|) and the definition of the critical oscillators, one has never 

T {N C ) > for ? - > /c(at c ).) It follows 
from flJJJ) and from T" < that V(£-j) < r'(tic(N c ),j) < 0. In addition, {TU implies 



Then, we obtain T^ Nc) - r'(^-j)(^ - ^j) < 0. 

• If the term -r"(£_ 3 -)(0j - %pj) > 0, it is associated with the term j ^ /C(l). (According to (Tlj) 

and the definition of the critical oscillators, one has never — r'(£_j)(0j- — V'j) > for j < /C(l).) It 
follows from JTSJ and from T" < that r'(£ K(1)j ) < r'(£_j) < 0. In addition, JJJ]) implies 

-(-I) 1 (^(i) ~ - fa - < ~ $j) < • 
Then, we obtain -T'^^fa - i>j) + rj 1} < 0. 

Finally, every term of ©-© is either negative or can be associated with a unique term so that the 
addition of both is negative. Thus, one has 

!ll®-*ll (1) = E r(fe) <o- 
fc=i 

For the other situations (r" > 0, V > 0), the proof follows on similar lines (except that the terms of 
might be associated with the terms of T( fc -i) instead). 

□ 



G 



Theorem 13.11 proves the contraction property inside the open cone C, but not for an orbit that lies in 
the boundary of C. However, the contraction property also holds for each of the (N — 2)-dimensional cones 
C[ = {& £ C\8i-i = 6i}, i = 1, . . . , N (with 9q = and On = 27r), that partition the boundary of C. A cone 
C[ corresponds to the synchronization of two oscillators and is therefore invariant under the dynamics ^ 
— since two synchronized oscillators remain synchronized forever. The result on the contraction of ([3J in 
the open cone C[ — defined as the interior of C[ — is summarized in the following corollary. 

Corollary 3.1. Under the same assumptions as Theorem \3.1l the dynamics (O defined in a cone C[ (i = 
1, . . . ,N) are either contracting (ifV <Q) or expanding (ifV > 0) in the open cone C[ with respect to ©. 

Proof. The proof is a straightforward corollary of the proof of Theorem 13.11 

In the cone C' x . Since 6\ — 0, the dynamics (O are reduced to the (N — 2)-dimcnsional dynamics 

JV-l 2V-1 

§ k = 2T(0 k ) + Y F0k - Oj) - Y r (-^i) - r (°) k = 2,...,N-l, 

3=2 3=2 
3¥fc 



which are similar to the dynamics in dimension N — 2. Then, the proof follows on similar lines as 

rjcL in is replaced by 2T$L 
r (k) 
-K(k) 



the proof of Theorem 13.11 The only difference is that the term TJ£L in is replaced by 2 TwL , but 



this modification docs not affect the validity of the proof. (In the proof, each term 2 is exclusively 

(fe+i) 



associated with another term 2T f ~^ k+1 yj 

In the cone C[, i = 2, . . . , N — 1. Since 0j = 6%-\, the dynamics Q arc reduced to the (N — 2)-dimensional 
dynamics 

N-l N-2 
3 = 1 3=1 

N-l N-2 

§, = r(e l )+ Y r (^-^) + r (o)- Y r(-§j)-2r(-§i) 

3 = 1 3=1 

j=£{i-l,i} 

and the proof follows on similar lines as the proof of Theorem 13. II The terms t[ in ([7]), with IC(k) ^ i, are 
replaced by 2 T^' and the term r"(£_i)(0j — ipi) ofT^ is replaced by 2r'(£_j)(0j — ipi). The reader can easily 
verify that these modifications do not affect the validity of the proof. (In the proof, the above-mentioned 
terms are exclusively combined between themselves.) 

In the cone C^. The situation is similar to the case of the cone C[. □ 

In addition, an equivalent result also holds for every intersection of several cones C' i: intersection that 
corresponds to the synchronization of three or more oscillators and that is therefore invariant under the 
dynamics . The result is summarized in the following corollary. (The proof follows on similar lines as the 
proof of Corollary 13.11 and is not detailed in the present paper.) 

Corollary 3.2. Let I denote a non-empty subset of {1, . . . , ./V}. Under the same assumptions as Theorem 
\3.1[ the dynamics §3§ defined in the intersection C'j = Hie/ C'i are either contracting (ifV < 0) or expanding 
(ifV > 0) in the interior of C'j with respect to @ . 

4. Collective behaviors of monotone oscillators 

Theorem 13.11 provides a contraction property for networks of monotone oscillators. This property is now 
shown to determine the asymptotic behavior of the network. 
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Theorem 4.1. Consider a network of N monotone phase-coupled oscillators with the dynamics ([T]). Pro- 
vided that the coupling function has a curvature of constant sign on (0, 2n), for almost every initial condition, 

• the oscillators asymptotically converge to the incoherent splay state configuration 0* k — k^r, with 
k G {1, . . . , N — 1}, if the coupling function is strictly decreasing on (0, 2tt); 

• the oscillators achieve perfect synchronization 0*(t) = 6%(t) = ■ ■ ■ = 0%(t) (i.e. 0% G {0, 27r} ) in finite 
time if the coupling function is strictly increasing on (0, 2n). 

Proof. Convergence to the splay state if V < 0. We first note that an orbit <fr(0, t) of with an 
initial condition in the open cone C cannot reach the boundary of C. Indeed, straightforward computations 
show the repelling property of the boundary: 

0i «r(0+) -r(27r") > for 0i ret), 

4 +1 -4«r(o+)-r(27r-)>o loth^Sk+i, (15) 

^at-i wr(27r-)-r(0+) <0 for (9jv-i ~ 2tt. 

Then, Theorem 13.11 implies that the dynamics ^ are contracting in C for all t and, for any small fixed At, 
the discrete-time mapping i— >• $(0, At) is also contracting in C. From the contraction mapping theorem, 
it follows that, for all initial conditions in C (the interior of the cone), the solutions of asymptotically 
converge to the unique fixed point lying in C. This unique fixed point corresponds to the splay state 0* k = k^S-, 
which concludes the first part of the proof. 

Convergence to synchronization when V > 0. Consider an orbit <&(0,t) of ([3]), with G C that is 
not the fixed point 0% = Theorem 13.11 implies that the distance between the orbit and the fixed point 
monotonically increases. As a consequence, the orbit converges toward the boundary of C (corresponding 
to the synchronization of at least two oscillators). In addition, (|15|) shows that the orbit approaches the 
boundary with the finite velocity |r(27r _ ) — T(0 + )| and therefore reaches the boundary in finite time. 

The orbit subsequently evolves within a (N — 2)-dimensional invariant cone C[. Given Corollary 13.11 
the reduced dynamics in C[ satisfies the contraction property. It follows that any fixed point lying in C[ is 
unstable, so that, for almost every initial condition, the orbit could not have reached C[ through a fixed 
point. Next, the same argument as above proves that the orbit reaches in finite time the boundary of C[, an 
event that corresponds to another finite-time synchronization of (at least) two oscillators. 

The orbit subsequently evolves within a (N — 3)-dimcnsional cone C\. Given Corollary 13.21 the reduced 
dynamics in C\ still satisfies the contraction property, so that the above argument can be repeated. 

Finally, the argument is repeated as many times as a pairwise synchronization occurs and after (at most) 
N — 1 successive pairwise synchronizations, the network achieves full synchronization in finite time. □ 

The two asymptotic collective behaviors of monotone oscillators are illustrated in Figure 0] Monotone 
phase-coupled oscillators obtained through the averaging of pulse-coupled leaky intcgrate-and-firc oscillators 
(see Q) satisfy the mild curvature assumption on T, so that the result applies. 

Theorem 14.11 is an illustration that contraction has strong implications on the asymptotic behavior. 
The global behavior of a network of monotone oscillators is in fact reminiscent of the global behavior of 
Kuramoto model, which is a gradient system. An important feature of the present model with respect to 
Kuramoto model is brought up by the discontinuity of the coupling function: the synchronization takes 
place in finite time, and the splay state is an isolated fixed point, whereas it is a N — 3 dimensional manifold 
in Kuramoto model. A consequence of that difference is that the asymptotic behavior of (pTJ) is robust to 
small heterogeneity in the natural frequencies u> in firing monotone oscillators Q whereas the asymptotic 
dynamics of Kuramoto model can be highly complex even for small heterogeneities [lij . 

5. Conclusion 

We have investigated the global stability properties of populations of monotone phase-coupled oscillators, 
thereby complementing local results presented in the earlier study In particular, we showed that mono- 
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°Q 20 40 ^30 ^ 5 ~l0 15 2fT 25 

Time Time 

(a) Incoherent splay state (b) Synchronization 

Figure 4: N = 10 monotone phase-coupled oscillators exhibit two opposite asymptotic behaviors, (a) With the decreasing 
coupling function r(9) = 1/7V(0.1 + e~ 9 ), the oscillators asymptotically converge to the splay state, (b) With the increas- 
ing coupling function T(6) = —1/N(0.1 + e~ e ), several pairwise synchronizations occur, until the network achieves perfect 
synchronization in finite time. 



tone oscillators only display two asymptotic collective behaviors: for almost all initial conditions, monotone 
oscillators either synchronize in finite time or asymptotically converge to the unique splay configuration. 

The global stability analysis relies on a strong contraction property of the dynamics, which is the main 
result of this paper. Interestingly, the contraction is not captured with respect to a quadratic norm, but 
through a 1-norm that has the interpretation of a total variation distance. In a general context, this result 
stresses the key role of 1-norms on cones to connect the monotonicity of a system to its stability properties. 
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